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We extend the SU(4) model [1-5] for high-T c superconductivity to an SU(4)£ model that permits explicit 
momentum (k) dependence in predicted observables. We derive and solve gap equations that depend on k, 
temperature, and doping from the SU(4)j. coherent states, and show that the new SU(4)j. model reduces to the 
original SU(4) model for observables that do not depend explicitly on momentum. The results of the SU(4)j. 
model are relevant for experiments such as ARPES that detect explicitly k-dependent properties. The present 
SU(4)t model describes quantitatively the pseudogap temperature scale and may explain why the ARPES- 
measured T* along the anti-nodal direction is larger than other measurements that do not resolve momentum. 
It also provides an immediate microscopic explanation for Fermi arcs observed in the pseudogap region. In 
addition, the model leads to a prediction that even in the underdoped regime, there exist doping-dependent 
windows around nodal points in the /c-space, where antiferromagnetism may be completely suppressed for all 
doping fractions, permitting pure superconducting states to exist. 

PACS numbers: 71.10.-w, 71.27.+a, 74.72.-h 



I. INTRODUCTION 

The mechanism that leads to high-temperature supercon- 
ductivity in the cuprates remains an open question despite 
intense study for the past two decades. Although the field 
has been challenged by many high-quality data from differ- 
ent types of measurement, there is no uniformly-accepted the- 
oretical picture that can offer a unified and consistent de- 
scription for these data. It is often believed that underlying 
physics can be understood in terms of individual particles 
interacting through appropriately-chosen interactions. How- 
ever, even with greatly-simplified Hamiltonians, describing 
collective motion in these strongly-correlated, many-electron 
systems has had only limited success. This has led some 
authors [6] to conclude that the many-body correlations in 
cuprates are so strong that dynamics may no longer be de- 
scribed meaningfully in terms of electrons and must be de- 
scribed instead in terms of new effective building blocks that 
fractionalize spin and charge. 

Simplicity is a kind of beauty in physics. Even if one 
could solve the problem with the help of large-scale numer- 
ical calculations, such a practice may not be interesting be- 
cause physics can be completely buried in the numerous con- 
figurations used in the calculation. On the other hand, al- 
ternative approaches to many-body problems have been pro- 
posed. One is the method of fermion dynamical symmetries 
[7]. This approach is based on the fact that collective motions 
in strongly correlated many-body systems are often governed 
by only a few collective degrees of freedom, and a quantum 
system exhibiting dynamical symmetries usually contains two 
or more competing collective modes. Once these degrees of 
freedom are identified and properly incorporated into a model, 
the problem may be considerably simplified and, most impor- 
tantly, the physics in such approaches may become transpar- 
ent. 

This is the philosophy that the SU(4) model of high- 
temperature superconductivity [1-5] is based on. For cuprate 



systems we have proposed that the most relevant collective 
degrees of freedom are af-wave superconductivity (SC) and 
antiferromagnetism (AF), and that coherent pairs (not individ- 
ual particles), formed from two electrons (or holes) centered 
on adjacent lattice sites, are appropriate dynamical building 
blocks of the wave function. The choice of this space, which 
is small in size but rich in physics, corresponds to a physically- 
motivated truncation of the huge Hilbert space corresponding 
to the original problem. 

It has been found [1] that these spin-singlet (D) and the 
spin-triplet (it) pair operators, when supplemented with the 
particle-hole type operators for staggered magnetization (Q), 
spin (S), and charge, constitute a 16-element operator set that 
is closed under a U(4) 3 U(l) x SU(4) algebra if the d- 
wave formfactor g(k) in D and % pair operators is replaced 
by sgn[g(Jb)]. The U(l) factor corresponds to a charge density 
wave that is independent of the SU(4) subspace because of the 
direct product. This implies that in the minimal U(4) model 
charge density waves do not influence the AF-SC competi- 
tion in lowest order and our discussions have been focused 
in the SU(4) subspace with its coherent-state approximation 
[1, 2], It has been further discovered [3] that the SU(4) sym- 
metry is a consequence of non-double-occupancy — the con- 
straint that each lattice site cannot have more than one va- 
lence electron. This suggests a fundamental relationship be- 
tween SU(4) symmetry and Mott-insulator normal states at 
half-filling for cuprate superconductors. 

Thus, the SU(4) model has ingredients of competing AF 
and SC modes, if-wave D and % pairs (entering as "preformed 
pairs" that are mixtures of the two kinds of paired states under 
the SU(4) constraint [4]) in wave functions, and non-double- 
occupancy imposed by the symmetry [3]. All of these appear 
to be relevant for the physics of cuprate superconductors. For 
data that do not resolve an explicit k dependence, the coherent- 
state solutions of the SU(4) model (with properly adjusted pa- 
rameters for effective interaction strengths) can consistently 
describe SC gaps, pseudogaps (PG), and the corresponding 
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transition temperatures T c and T* in cuprates, as demonstrated 
inRef. [5]. 

However, there are experimental indications for explicit k- 
dependence that are observed by experiments such as angle- 
resolved photoemission spectroscopy (ARPES) [8, 9]. These 
data probe electrons near the Fermi surface having particular 
k directions. To describe ^-dependence in energy gaps, we 
must extend our original SU(4) model by displaying explicit 
^-dependence in the gap equations and their corresponding so- 
lutions. This is the goal of the present paper. 

The paper is organized as follows. In Sec. II, we out- 
line the SU(4) background by pointing out the assumptions 
made when the original ^-independent SU(4) model was con- 
structed. Sections III and IV are respectively devoted to 
presentation of the new fc-dependent SU(4),t model and the 
fe-dependent gap equations obtained using the generalized 
coherent-state method. We solve these gap equations in Sec. V 
and give analytical solutions for the superconducting gap and 
the pseudogap. Finally, we discuss some immediate conse- 
quences of the SU(4),t model in Sec. VI, and a short summary 
is given in Sec. VII. 

II. DYNAMICAL SYMMETRIES AND THE ORIGINAL 
SU(4) MODEL 

Interactions in dynamical symmetry theories are deter- 
mined by symmetry groups [1]. A general SU(4) Hamiltonian 
with pairing and AF interactions can be written as [10] 

H = H -V d -V 7[ -V q , (1) 

where Hq is the single particle (s.p.) energy, and V d , V n , and 
V q are the two-body spin-singlet pairing, spin-triplet pairing, 
and AF interactions, respectively: 

Ho = Y. £ k n k (2a) 

k 

V * = Y,Gl k ,D\k)D(k') (2b) 
k,k' 

v * = L G **'* t (*)-*(*') (2c) 

k,kl 

v g = l>^e(*)-e(*')- (2d) 

k,k! 

The operators appearing in Eqs. (2) can be expressed as 

D\k) = g(k)cl/_ kl (3a) 
= g{k)cl +q .c\. (3b) 

Gy(*) = 4 +9 /*J> ( 3c > 

where n}:(k) and Qij{k) are, respectively, tensor forms of 
7t^(k) and Q(k). In Eqs. (3), . creates an electron of mo- 
mentum k and spin projection i,j = 1 or 2 (= | or J.), and 
q = (n, K, n) is an AF ordering vector. The d-wave formfac- 
tor, 

g(k)=g(k x ,ky) = cosk x -cosk y , (4) 



appears in (3a) and (3b) because of strong experimental evi- 
dence that in cuprates the coherent pairs exhibit d-wave orbital 
symmetry [11]. Energy gaps thus generally are fc-dependent 

M*) = L G °kA D '( k )) ( 5a ) 
k' 

M*) = L^, «(*')) (5b) 
k' 

A #) = L4k'(Q^ k '))- ( 5c > 
ki 

The discussion to this point is general and no approxima- 
tions have been made. In our original fe-independent SU(4) 
model [1-5], we have introduced approximations through the 
following assumptions 

g(k) » sgn[g(k)} (6a) 
£k ~ £ (6b) 
G' kk , w G (i = 0,l), X° kk '~X°- (6c) 

Assumption (6a) removes the ^-dependence from formfactors 
in the pair operators, and assumptions (6b) and (6c), respec- 
tively, replace the s.p. energy and interaction strengths with 
fc-independent constants. These approximations thus lead to 
fc-independent gaps 

A d = (?(rf) A n = G l (nl) A q = X°(Qz), (7) 

which are expressed in terms of the collective operators 

D + =£sgn[g(*)]4/_ ti 

k 

4 = L s ^[g(k)]cl +q / kj 

k 

k 

The preceding equations constitute the basis for discussions 
in Refs. [1-5], and all our previous SU(4) results are obtained 
within this framework. As most cuprate data presumably 
represent weighted averages over contributions of different k 
components, the original SU(4) scheme works well. In Ref. 
[4], we derived and solved fe-independent (but temperature 
and hole-doping dependent) SU(4) gap equations, and used 
the results to construct generic gap and phase diagrams. We 
compared the results with some representative cuprate data in 
Ref. [5] and found that, for data that do not resolve an ex- 
plicit k dependence, the coherent-state solutions of the orig- 
inal SU(4) model can consistently describe SC gaps, pseu- 
dogaps, and the corresponding transition temperatures T c and 
T*. 



III. THE ^-DEPENDENT SU(4)* MODEL 

As we have noted, there is experimental evidence for ex- 
plicit fc-dependence of certain observables in the cuprates. Al- 
though interpretation of some of the results remains somewhat 
controversial, their momentum-dependent nature is clear. One 
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FIG. 1 : Geometry and definitions in the k x -k Y plane, where k is hole 
momentum. 



gok = 1 - cos k 




FIG. 2: The curve \g(k)\ and its maximum value goi = 1 ~ cosk for 
a given momentum k = (k x ,k y ), under the constraint (11). 

example is the observation of Fermi arcs in angle-resolved 
photoemission spectroscopy (ARPES) data [8, 9]: ARPES 
measurements suggest that the Fermi surface is gapped out 
arc wise in the pseudogap region below T*, indicating clear 
anisotropy of the PG in the A:-space. 

In order to describe momentum dependence of energy gaps, 
we must extend our original SU(4) model [1-5] in a way that 
restores the ^-dependence that is washed out by the assump- 
tions in Eqs. (6), but preserves the SU(4) symmetry. The re- 
placement (6a) for the pair operators is a necessary condition 
for preserving the SU(4) algebra [1], which is required physi- 
cally because it imposes the non-double-occupancy condition 
[3]. Therefore, the only way to restore ^-dependence in en- 
ergy gaps but keep the SU(4) symmetry (and its associated 
non-double-occupancy constraint) is to modify (6b) and (6c) 
to allow the s.p. energy £ and the interaction strengths to carry 
^-dependence. The s.p. energy term e is less important in this 
regard because it does not contribute to energy gaps and tran- 
sition temperature in our formalism [4] . We may thus employ 
it in the most general form 

Without loss of generality, g(k) can be written as the prod- 
uct of the absolute value and a sign 

g(k) = \g(k)\x S gn[g(k)}. (9) 

Therefore, the approximation (6a) implies that in our original 
SU(4) model we have assumed that the magnitude |g(Jfc)| is 



unity, regardless of k. (Note that this is also the condition used 
to close the algebra of the SO(5) model [12, 13].) Instead of 
\g(k) \ = 1, we now introduce 

\g{k)\na k = g 0k 8(e k ), (10) 

where go<r is the maximum value of |g(Jfc)|. In Figs. 1 and 
2, we illustrate the associated geometry and definitions. In 
our notation, k = (k x , k y ) is the electron momentum under the 
constraint 

{7t-k x ) 2 + {7t-k y ) 2 =k 2 , (11) 

where k is the hole momentum with k its azimuthal angle, as 
shown in Fig. 1. In Eq. (10), 8(0^) takes a value of unity ex- 
cept in a narrow region around the nodal points (correspond- 
ing to \k x \ = \k y \ or Ok = ±n/4 for the first Brillouin zone; see 
Fig. 1), where it quickly diminishes and vanishes exactly at 
the nodal points. A possible mathematical expression could 
be of the Gaussian type 

( ?!C*1± \ 2 

8(6 k ) = 1 - e \ _5B 7 with A0 < 

where AO measures width of the Gaussian. With very small 
AO, the exponential term has a negligible contribution to the 
average, which ensures the averaged 8(0^) equal to 1. The so- 
defined 8 (Ok) becomes exactly zero at 0^ = ±7r/4. Therefore, 
our pairing gaps have nodes at Ok = ±7t/4, which agrees with 
experiments. 

The behavior of \g(k) \ is illustrated in Fig. 2. It is easy to 
show that 

g0k= |l — cosfe| . (12) 

Thus, for the first Brillouin zone k x and k y can take values 
from zero to n, while go*.- changes from to 2. (The assump- 
tion \g(k) \ — 1 in the original SU(4) model is thus equivalent 
to taking an average of go* over k.) Equation (10), with its ex- 
plicit dependence on k, clearly improves on the original SU(4) 
model for observables having a possible k dependence. 

With the new approximation (10) for the pairing in- 

teraction strengths in Eq. (2) are 

G^ = Gfoi^ (» = 0,1). (13) 

The factor sgn[g(Jfc)] in the product (9) remains unchanged in 
the pair operators, which ensures preservation of the SU(4) 
symmetry. 

The ^-dependence of the AF interaction % kk i follows from 
the nature of exchange interactions. The corresponding ma- 
trix elements are proportional to the wavefunction overlap 
between the states, which are one-particle, one-hole states 
with momenta (k + q,k). The of-wave symmetry in the pair 
structure implies that the amplitude of a pair wavefunction 
with two electrons having momenta (k, —k) in the background 
mean-field of the SU(4) collective subspace is proportional to 
g(k). This means physically that the two electrons in a pair 
favor aligning their momenta k along the Cu-O bond direc- 
tion (maximum of |g(Jfc)|) rather than along the diagonal to 
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the Cu-O bonds (nodal direction of g(k); see Fig. 1). Since 
the wavefunction for momentum fe + q differs from that for fe 
only by a sign [1], and a hole wavefunction is the conjugate 
of a particle wavefunction, the amplitude of a particle-hole 
wavefunction with momenta (Jfc + ^,Jt) should be very similar 
to that of a pair and thus proportional to g(k). Therefore, we 
can write 



Xu? = x g kg a, 



(14) 



where g k = \g(k)\. Note that for the g k factor in Eq. (14), no 
approximation like the one in Eq. (10) is necessary. Explicitly, 
we mean here that g k = \g(k) \ = \ cosk x — cos fey |. 

Inserting Eqs. (13) and (14) into (2), we can rewrite the 
Hamiltonian (1) as 



k k.k 1 

-G\^a k a k ,KHk)-K{k') 

k,k' 

-X°L8k8k'Q(k)-Q(k'). 
k,k' 



(15) 



The fe-dependent Hamiltonian (15) possesses a ^0 SU(4)j 
symmetry, with 15 fe-dependent generators 

D + (fe) = D + (fe)+D + (-fe) 

* + (Jfc) = ^ t (fe) + ?r t (-fe) 

Q(k) = Q(k) + Q(-k) 

M(jfc) = Af(Jt)+Af(-Jt) 

S(k) = S(k) + S(-k), 

where M(fe) and S(k) are, respectively, the charge and the spin 
operators. For each k, the commutation relation among gener- 
ators, the structure of subgroup chains, and their correspond- 
ing properties are analogous to those of the original SU(4) 
group structure [1]. We term this fe-dependent extension [14] 
of the original SU(4) model the SU(4)* model. 

IV. ^-DEPENDENT GAP EQUATIONS 

In the preceding section, we demonstrated that with a better 
approximation to the absolute value of the formfactorg(fe) it is 
possible to introduce explicit fe dependence through a symme- 
try structure that corresponds to a product of SU(4) groups, 
each labeled by fe. Therefore, the following discussions for 
gap equations and their solutions for a given fe are rather sim- 
ilar to those in Ref. [4] for the fe-independent SU(4) model. 

By analogy with discussions in Appendix B of Ref. [4], 
under the coherent-state (symmetry-constrained, generalized 
Hartree-Fock-Bogoliubov) approximation, one obtains for 
the fe-dependent case 



2k*± Vfct (fifct - A ) - Afct {u\ ± - 



■vL) 







(16) 



with 



and 

A d (k) = Gg<7 k £^(Dt(fc')) 

k'>0 

A n {k) = G?a*E^ (*!(*')) 

k'>0 

\(k) = zVE^(Qz(*'))- 

k!>o 

fe' > in the above and following equations means k' x > or 
k' y > 0. Solving Eq. (16) gives fe-dependent occupation proba- 
bilities 



2 

u k± 



1 



ek± 



i - 



<?k± 



and a quasiparticle energy 

e k ± = yW-A) 2 + A2 ± . 
The gap equations in fe-space can then be obtained: 
G° Q a k 



A d (k) 
A K (k) 
\(k) 

-2x 



E a k'( w k'+ A k'++ w k'- A k'- 

f>0 



(17a) 



G\a k 



X°gk 



f>0 



V k'+ A k'+- W k'- A k'-) (l 7b ) 



E^K' + [A,(fe')+^] 



k'>0 

+w k ,_[A q {k')-^W 

|EK' + K( fc, )+^] 



(17c) 



(17d) 



with 



w k ± 



Pk±(T) 



e k ± 

P k± (T)=tmh 



K — & - £ k 

e k ± 
2k B T 



£k± = £k T A q (fe) A k± = A d (fe) ± A * (fe) 



In Eq. (17d), Q. = Y, k >o i s tne maximum number of doped 
holes (or doped electrons for electron-doped compounds) that 
can form coherent pairs, assuming the normal state (at half 
filling) to be the vacuum, x is the relative doping fraction in 
the model [4]. Positive x represents the case of hole doping, 
with x = corresponding to half filling (no doping) and x = 1 
to maximal hole doping. The true doping P is related to x by 
jc~4P[4]. 

The fe-dependent gap equations (17) are coupled algebraic 
equations. By solving these equations, one can in principle 
obtain fe-dependent (also temperature and hole-doping depen- 
dent) energy gaps. However, general and exact solutions are 
difficult because gaps for given fe are related to all other fe 
points, which means that solutions at each fe point are not in- 
dependent from the other fe-components. In the next section, 
we show that one can obtain analytical solutions by applying 
some approximations. 
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V. SOLUTIONS FOR jfc-DEPENDENT GAP EQUATIONS 



which implies that 



We may greatly simplify the solution of Eqs. (17) through 
the following three steps. First, we replace the quantities in 
the summations on the right hand side of (17) with their cor- 
responding mean values: 



Ajtt => A t (t = +,-,q) 

with 

P±(T) = tanh 



XL 



w± = — , 



( e ± 



\2k B T 



e± = V(A 7 ±A,) 2 - 



-r2 



The functions a k and g k in the summations can then be sim- 
plified as 



E °k< 

*'>o 

E sa 

k'>0 



780 



E so 

Jf>0 

L 8=^-8- 
k'>o L 



(18) 
(19) 



The second level of simplification is based on physical con- 
siderations. Experimentally-measured energy gaps are dom- 
inated by contributions from near the Fermi surface. There- 
fore, we assume that measured gaps may be approximated by 
their values at k = kf. Using this approximation and the aver- 
age values introduced in the first approximation step, we can 
write for the gap equations of (17) evaluated at k = kf. 



A d (k) = 


Q 
4~ 


A K (k) = 


Q 


4" 




Q 


A q (k) = 


"4 



-2x = 



^G°i8o8oS(9k f ) (w+A+ - w 

Ao?7(^){w+[A 9 + A 7 ] 
+W_[A,-X 7 ]} 

W+fi q +17\-W-\A q -T'\, 



where g Q = g Qkp and 



7(0%) 



8{k) 



80 



(20a) 
(20b) 



(20c) 
(20d) 



(21) 



Equations (20) are ^-dependent gap equations constrained on 
the Fermi surface through 



{n - k x f + (tz - k y 



(22) 



It can be shown that y(9 kf ) is independent of \kf\, to high ac- 
curacy, and therefore can be considered in later discussions to 
be a function of azimuthal angle only. 

In the third simplification step, we assume the average val- 
ues A T and A' to be proportional to the unknown quantities 
A T (Jt) and X k , respectively, with a constant of proportionality 
R: 



A T = RA T (k) (t: 



,q) l' = Rli (23) 



e±=Re k± . 



(24) 



The parameter R serves as a renormalization factor that cor- 
rects on average for the errors caused by the approximation 
and is determined by fitting data. With (23), Eqs. (20) now 
become 



-2x 



-008(6^) [w+A+(*)+w_A_(*)] (25a) 
jGi8(6^)[w + A+(k)-w-A-(k)] (25b) 
{w+[A q (k)+X k ] 



= TX 



4 '" 7 

+w-[A q (k)-^}} 
w + [A q (k)+X' k ] 



(25c) 
(25d) 



with 



and 



Gi = GUogo 



X = X°8oY8 



w± = 



P±(T) = tanh 



Re k± 
2k B T 



(26) 



In the above equations, 7 is the average value of 7(0^). 

The simplified gap equations (25) can now be solved ana- 
lytically. They have the same structure as the gap equations 
discussed in the ^-independent SU(4) model [4], except that 
the interaction strengths in the present case are ^-anisotropic. 
Therefore, all the SU(4) formulas in Sections III- VI of Ref. 
[4] remain valid, provided that the following replacements are 
made for the singlet-pairing, triplet-pairing, and antiferromag- 
netic coupling strengths, respectively: 

Go^Go8(9 kf ) Gi-GiSftfcf) X^XY(0 k{ )/7- (27) 

For example, if we introduce the the doping parameter x de- 
fined in §II.c of Ref. [4], the fe-dependent critical hole-doping 
fraction is [compare Eq. (23) of Ref. [4]] 



Xqd = x q(9kf) 



IXY(0k f )/Y-G Q 8(e k{ ) 

xrio^/r-G^y 



(28) 



and the T — energy gaps at the Fermi momentum kf for x < 
x q Q are obtained as 



V*) 

K 



Go8{9 kf )^x(xJ -x) 

LI I 

-^G ] _8(d kf )^x{x qe -x) 



(29a) 
(29b) 
(29c) 



~2 



X^-G^) 



XqB (l ~X q ex) 



-—G]_8(e kf )x, 



(29d) 
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while for x > x q g we obtain A q (k) = A K (k) = and 



K 



(30a) 
(30b) 



for the solutions. The energy gaps obtained in Eqs. (29)-(30) 
are ^-anisotropic. The pairing gaps have nodal points at k x — 
k y , where S(0 kf ) = 0. The pseudogap A q (k) is a function of k 
by virtue of the factor y(9 kf ) defined in Eq. (21). 

Because all SU(4) formulas in Sections III-VI of Ref. [4] 
remain valid, it is easily proven that the PG closure tempera- 
ture T* acquires the same g{k) dependence as the pseudogap, 
and we obtain for the PG closure temperature 



r ( *)=*^* (1 -* 2) 



4fe R 



(31) 



We do not expect a corresponding effect in the superconduct- 
ing region because below T c the pairing gap opens and the en- 
tire Fermi surface will be destroyed except at the nodal points. 
We find a superconducting transition temperature 



(32) 



which has no g{k) factor. 

VI. DISCUSSIONS AND PREDICTIONS 



Most experimental techniques do not resolve k and we ex- 
pect for those that transition temperatures are dominated by 
contributions from near the Fermi surface (k = kf), averaged 
over all fe-directions. If one takes the average over 9 kf , then 

5(^)-i y(e kf )/y^i, 

and the gap equations of the ^-dependent SU(4)^ model and 
their solutions become identical to those obtained for the orig- 
inal ^-independent SU(4) model [4]. Thus our original SU(4) 
model predicts [4, 5] values of energy gaps and the corre- 
sponding transition temperatures that are (perhaps weighted) 
averages over k. These are relevant for comparison with ex- 
periments that do not resolve k. However, the explicit appear- 
ance of the anisotropic factor y(0 kf )/y in the gap solutions 
of the SU(4)£ model leads to some interesting new conse- 
quences. We note that although our following discussions are 
made through the anisotropic factor y(9 kf )/y, its relation with 
8(0 kf ) guarantees that the model still preserves the <f-wave 
nature and has nodes in the pairing gaps. In this section we 
discuss three predictions following from the new formalism 
that could have important implications for experiments that 
detect explicitly fe-dependent properties. 



A. Two pseudogap closure temperatures: the maximum and 
the averaged 

The ^-dependent PG closure temperature T*(k) in Eq. (31) 
differs from the A:-averaged one derived in Eq. (49) of Ref. [4], 



4k B 



(33) 



by the factor y{9 kf )/y. We know that y{9 k{ ), and also T*(k), 
take their maximum values at the antinodal points; for exam- 
ple, 

7max(^ f )|e % =0,f = 1 

(see Figs. 1 and 2). We can denote the maximum PG closure 
temperature as 7^ ax . Thus, the SU(4)^ model predicts two PG 
closure temperatures that are related to each other through 



T* = T* ly 

max -'ave/ /■ 



(34) 



It is straightforward to evaluate the averaged y value by inte- 
gration. Noting that Eqs. (29) are restricted to x < x qd , we can 
define the maximum allowed azimuthal angle C through the 
condition x = x q g(0 c ). We then have 



2 f(> 
'■= - / 
n Jo 



8(k(0)) 



go 



dd. 



(35) 



In the above calculation, we have used for the integrand the 
expression (21) for y(9 kf ) and the constraint (22). The result- 
ing y depends on the size of the Fermi surface \kf\. Assuming 
an isotropic hole Fermi surface, we have 

kt 2 =2n{l+P). 

Therefore, y is essentially a hole-doping-dependent quantity. 
In Fig. 3, we show the behavior of 1 jy as a function of doping 
P assuming coupling-strength parameters characteristic of the 
cuprate superconductors. As one can see, it has a nonlinear 
dependence on doping, taking the maximum value 1 .6 at very 
small dopings, falling rapidly between P = 0.05 and 0.08, and 
continuously decreases but with a smaller rate until it reaches 
unity at the critical doping P = 0.18. 

We thus obtain two distinct PG closure temperatures, r a * e 
and r^ ax , having the same microscopic origin [4, 5], but dif- 
fering in the kinds of experimental observables for which 
they are appropriate. The largest difference between the two 
is found for small doping; they take similar values at large 
dopings, becoming identical at the optimal doping point. In 
Ref. [5], experimental values of T c and T* [15] that do not re- 
solve k were compared to our theoretical T c and T* v e . In Fig. 4, 
we re-plot these values in green for T c (with open triangles for 
data and dotted curve for theory) and blue for T* vss (with open 
squares for data and solid curve for theory). Above the blue 
(solid) curve, we now add the maximum PG closure tempera- 
tures r^ ax in red (with open circles for data and dashed curve 
for theory). Because of the 1 jy factor shown in Fig. 3, the 
red (dashed) curve lies well above the blue (solid) curve at 
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FIG. 3: The doping-dependent 1/y factor. The calculation employs 
Eq. (35) and utilizes realistic interaction strengths x, Go and G\ taken 
from Ref. [5], with the pairing onset at P = 0.05. 
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FIG. 4: (Color online) SU(4) cuprate phase diagram compared with 
data. Strengths of the AF and singlet pairing correlations were deter- 
mined in Ref. [5] by global fits to cuprate data. The PG temperature 
is T* and the SC transition temperature is T c . The AF correlations 
vanish, leaving a pure singlet rf-wave condensate, above the criti- 
cal doping P q . Dominant correlations in each region are indicted by 
italic labels. Data in green (open triangles) and blue (open squares) 
are taken from Ref. [15], and those in red (open circles) from Ref. 
[16] (arrows indicate that the point is a lower limit). 



low dopings. In Ref. [16], Campuzano et al. reported their 
ARPES data (plotted as red circles in Fig. 4). In the un- 
derdoped regime it is clear that the results from Refs. [16] 
and [15] differ substantially. Because the ARPES experiment 
typically detects ^-dependent properties, we suggest that the 
data from Ref. [16] actually measure 7^ ax , as predicted in the 
present paper, while the data cited in Ref. [15] measure the k- 
averaged T* ve , as described in our earlier paper (Ref. [5]). We 
emphasize that in this interpretation the two types of experi- 
ments are seeing the same underlying physics, but the obser- 
vations differ because what is actually being measured differs 
in the two cases. 



B. Temperature-dependent Fermi arcs 

The pseudogap closure temperature T*(k) is anisotropic in 
k. Combining Eqs. (31), (33) and (34), we have 

T*{k) = T^y{e kf ) = ^\g{k i )l (36) 
go 

where the doping-dependent quantity Z^ ax is the maximum 
value of T*(k) in the antinodal direction, dj^ = or n/2. For 
an arbitrary temperature T < 7^ ax , the ^-dependent pseudogap 
closes when T = T*(k), which is equivalent to the require- 
ment that 

| cosfe v - cosfevl = g Q (T/T^ x ), (37) 

upon substituting (4) and (36). This equation says that the 
magnitude of the li-wave formfactor (4) that expresses the 
nodal structure [11] in cuprate superconductors is propor- 
tional to the scaled quantity T /7^ ax , with a proportionality 
factor go that is related to the size of the Fermi surface. 

Simultaneous solution of Eqs. (22) and (37) gives values of 
k x and k y where the pseudogap closes at a given P and T. Fig- 
ure 5 illustrates the solution of (22) and (37) graphically for 
several temperatures at fixed doping. The solution of Eq. (37) 
is represented by the curves bounding the shaded regions and 
the solution of Eq. (22) is represented by the Fermi surface 
curves in each corner. The intersection of these curves defines 
two simultaneous solutions in each of the four quadrants that 
bound the surviving part of the Fermi surface (heavier portions 
of the curves in Fig. 5). In the shaded portions of Fig. 5 the 
Fermi surface has been destroyed by the pseudogap, leaving 
only a vestigial Fermi arc between the shaded regions. 

The solution in Fig. 5 represents a derivation of Fermi-arc 
structure expected in the underdoped regime above T c . Below 
T c , the Fermi surface is completely destroyed by the opening 
of the pairing gap except at the nodal point, since the pairing 
gap has no y(0jt f ) dependence [see Eq. (29)]. We emphasize 
that the calculations presented in Fig. 5 do not involve any 
new parameters as long as the PG closure temperatures have 
been calculated (presented in Fig. 4). 

Kanigel et al. [9] have reported from their ARPES exper- 
iment arc lengths for slightly underdoped Bi2212. A direct 
reading of the Fermi arc length from Fig. 5 permits us to com- 
pare it quantitatively with the Kanigel data. Details will be 
published elsewhere [17]. 



C. Complete suppression of antiferromagnetism: pure 
superconducting states in underdoped compounds 

As we have discussed extensively in Ref. [4], the antifer- 
romagnetic correlation that plays a key role in understanding 
underdoped cuprates is completely suppressed at and beyond 
the critical doping point x q . A pure (af-wave) BCS supercon- 
ducting state occurs at zero temperature in the overdoped por- 
tion of the phase diagram. We now show that a similar situ- 
ation can also occur in certain fc-windows in the underdoped 
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FIG. 5: (Color online) Construction of Fermi arcs for doping P = 0.15 and values of T /T^x decreasing left to right. The ranges of k x and 
k y are — 71 to 71 and dashed lines indicate nodes. Hole Fermi surfaces in the absence of gaps are illustrated by full solid arcs in each corner. 
For T > T^ x a full Fermi surface exists; for T < 7^ ax , opening of the pseudogap destroys Fermi surfaces in the shaded regions (dotted lines), 
leaving arcs (solid lines) centered on the nodal lines. These arcs have absolute lengths that depend on P and T, but relative lengths that depend 
only very weakly on P and are determined almost entirely by the ratio r/r,* ax . The sizes of the shaded regions grow with decreasing T, so at 
very low temperature almost all of the Fermi surface becomes gapped and the Fermi arcs shrink to the nodal points as T /T^. dx — » 0. 
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quential figures, plotted for four different dopings in Fig. 7, 
show the behavior of the energy gaps as functions of the mo- 
mentum direction. Whenever 
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FIG. 6: The anisotropic factor y(6). In this figure, 6 Cl = 6 C , C , = 
71 / 2 — 9 C , and 6 C is determined by x c/ g(0 c ) = x. The value of y{8) 
at the critical angle C is denoted by y(0 c ) = yc(x), where y c (x) is a 
monotonically increasing function of doping x that becomes equal to 
unity when x = x q . 



regime, in which AF correlation is completely suppressed and 
a pure SC state emerges. 

This is another interesting consequence of the SU(4)j 
model due to the anisotropic factor y(fy f )/y. The critical 
doping point defined in Eq. (28), which is constant in the k- 
averaged SU(4) model [4], is now a function of momentum di- 
rection 8]^ because of the anisotropic factor y(6jt f )/y Conse- 
quently, for each given doping x there always exists a window 
in the momentum azimuthal angle, 9 C < 9k f < (%/2 — 9 c ), and 
centering at the nodal point 71 / 4, within which the AF correla- 
tion vanishes and only the pairing gap A c / exists. This follows 
because inside the window x q g < x; therefore, the solution 
(29) is not permitted but the solution (30) is. The critical an- 
gle 9 C is determined by the condition x q g(9 c ) = x. Figure 6 
illustrates the situation. 

Because y(% f )/y, and thus x q g, is a doping-dependent 
quantity, the above phenomenon depends on doping. The se- 
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9 C = 0, which means that there is no momentum space avail- 
able to Aq and A n , and the AF correlations and triplet pairing 
states are completely suppressed. Therefore, for x > x q the 
system can only be in a pure superconducting state at zero 
temperature. It can be seen that larger doping x implies a 
smaller critical angle C , and thus a wider pairing window, and 
that the width of the pairing window decreases rapidly toward 
zero as the doping goes to zero. 

In the original SU(4) model, we found that the critical dop- 
ing point defines a natural boundary (quantum phase transi- 
tion) between underdoped and overdoped regimes that have 
qualitatively different wavefunctions [4] . We termed the un- 
derdoped superconducting regime the AF+SC phase (anti- 
ferromagnetic superconducting phase); it is characterized by 
having all gaps nonzero but is dominated by AF and SC gaps. 
The present extension to the SU(4)^ model reveals the addi- 
tional feature that in this AF+SC phase the gaps are highly 
anisotropic in the momentum space, implying the possibility 
of a pure SC window around the nodal points even in the un- 
derdoped regime. The proposed existence of these pure su- 
perconducting windows may have considerable implication 
for the nature of the Fermi surface at low doping, for the 
Nernst effect, and for the relationship of impurities to inho- 
mogeneities in the underdoped region. These deserve further 
investigation. 

The above analysis has assumed T = for simplicity, but 
the basic picture should be valid also in the case with nonzero 
temperature. The formulation and solution of the gap equa- 
tions described in this paper can be extended to finite tem- 
perature using the methods described in Ref. [5]. While of 
considerable practical importance, this extension does not in- 
volve conceptually new ideas and will be deferred to a later 
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FIG. 7: (Color online) Dependence of the energy gaps Aj, Ajj, and A„ on the momentum direction 6 for several representative dopings P. For 
each doping there is a 0-window (yellow shaded region), centered at the node 6 = k/4 (indicated by the open circle), in which A 9 = A n — 0, 
and A ( / ^ 0. This window encompasses only a few percent of the Brillouin zone at low doping (tending to zero at zero doping), but rapidly 
expands to fill the entire Brillouin zone near the critical doping P ~ 0. 18 and beyond. 



paper. 

VII. SUMMARY 

In this paper, we have extended the SU(4) model for high-T^ 
superconductivity to include explicit momentum dependence 
in observables. To do so, we have started from a general 
SU(4) Hamiltonian and introduced a new approximation for 
the J-wave formfactor in the pair operators. This leads to 
the new SU(4)^ model, which retains explicit ^-dependence 
while preserving SU(4) symmetry. We have solved the gap 
equations derived from the SU(4)j coherent states with some 
plausible approximations, obtaining analytical solutions for k- 
dependent superconducting gaps, pseudogaps, and their tran- 
sition temperatures T c and T*. The new SU(4)j model re- 
duces to the original SU(4) model for observables that are 
averaged over all possible k directions. Therefore we pro- 
pose that the original SU(4) model describes the averaged fea- 
tures and thermal properties of cuprates, while the new SU^)^ 
model presented in this paper extends this description to de- 
tailed anisotropic properties in the fc-space. The present re- 
sults have been obtained for zero temperature but the formal- 
ism presented here may be extended to finite temperature in 
a manner similar to that extension of the ^-averaged SU(4) 
model. 

Because of an anisotropic factor y(0jt f )/y m the analytical 
gap solutions, the cuprate phase structure in the underdoped 
regime becomes even richer than that for fc-averaged observa- 
tions. We have discussed three immediate consequences that 
emerge in the new SU(4)j model: 

1 . We have suggested the possibility of two distinct, mea- 
surable, pseudogap closure temperatures: the maximum 
and the averaged. In the coherent-state SU(4) theory the 
pseudogap could be interpreted either as arising from 



competing AF and SC degrees of freedom, or alterna- 
tively as fluctuations of pairing subject to SU(4) con- 
straints [4, 5]. The proposed T* ve and 7^ ax share this 
same microscopic origin, but differ from each other by 
a doping-dependent factor. The temperature T a * e repre- 
sents PG closure temperatures that are averages over k, 
while r^ ax , which is generally higher than r a * e , repre- 
sents the pseudogap temperature expected if one retains 
explicit ^-dependence. Experimentally, then, we pre- 
dict that T* ve is the pseudogap temperature that should 
be measured in experiments that do not resolve k ex- 
plicitly, but (the generally higher) is the expected 
measured pseudogap temperature for experiments like 
ARPES that resolve k. 

2. We have provided a theoretical framework to under- 
stand ARPES Fermi-arc data. Using two analytical 
equations, we have obtained solution for the T /T* de- 
pendent Fermi-arc lengths in quantitative agreement 
with existing measurements. The essence of this result 
is the appearance of the new factor 7(0) in the PG clo- 
sure temperature that has been derived in Eq. (31). 

3. We have predicted the existence of doping-dependent 
windows in the momentum space where antiferromag- 
netic correlation is completely suppressed in the under- 
doped regime. Without AF competition, it is possible 
for pure superconducting states to emerge in these win- 
dows. Thus, we find that pure BCS-type superconduct- 
ing states can exist, not only in conventional supercon- 
ductors or in overdoped cuprate high-T c superconduc- 
tors (where such behavior is well established), but in 
localized islands even in underdoped cuprate supercon- 
ductors. It is of interest whether this prediction is re- 
lated to the recent observation of small pockets of well- 
defined fermi surface in underdoped cuprate supercon- 
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ductors. 

Some of these predictions (for example, the two temperature 
scales for pseudogap behavior) have the potential to recon- 
cile apparent discrepancies in existing data. All make predic- 
tions that can be tested in experiments capable of resolving 
^-dependent behavior. 

Finally, we note that the recent discovery of superconduc- 
tivity in layered iron-based transition metal oxypnictides [18] 
has generated a new wave of research interest. In place of 
copper and oxygen, the new compounds contain iron and ar- 
senic, and the highest critical temperature for them has already 
reached 55 kelvin [19]. It has been demonstrated in neutron- 



scattering experiments [20] that, like high-r t copper oxides, 
superconductivity in these iron-based materials is likely com- 
peting strongly with antiferromagnetic degrees of freedom. It 
will be of considerable interest to see whether approaches like 
the one presented in this paper, or other models capable of 
handling multiple competing degrees of freedom in strongly- 
correlated systems on an equal footing, can explain these new 
high temperature superconductors and their relationship to the 
old ones [21]. In particular, we note that what is already 
known about the new iron-based superconductors suggests 
that fc-dependent phenomena of the sort described in this paper 
should also be observable in these new superconductors. 
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